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1 Introduction 



One of the strongest constraints that a unified theory of the fundamental interactions 
must satisfy is that it leads to realistic cosmological models. There are known to 
be five consistent perturbative superstring theories in ten dimensions. (See, e.g., 
Ref. [|I|). The type II and heterotic theories each contain a Neveu-Schwarz/Neveu- 
Schwarz (NS-NS) sector of bosonic massless excitations that includes a scalar dilaton 
field, a graviton and an antisymmetric two-form potential. The interactions between 
these fields lead to significant deviations from the standard, hot big bang model based 
on conventional Einstein gravity and a study of the cosmological consequences of 
superstring theory is therefore important. 

In string cosmology, the dynamics of the universe below the string scale is deter- 
mined by the effective supergravity actions. The NS-NS string cosmologies contain 
the non-trivial fields discussed above. Recently (hereafter referred to as paper I), 
a complete qualitative analysis for the spatially fiat Friedmann-Robertson-Walker 
(FRW) and axisymmetric Bianchi type I NS-NS string cosmologies was presented. A 
central charge deficit was also included and found to have significant effects on the na- 
ture of the equihbrium points. This study unified and extended previous qualitative 
analyses of this system 

In Ref. (hereafter referred to as paper II), a phenomenological cosmological 
constant was introduced into the string frame effective action in such a way that it 
was not coupled directly to the dilaton field. Such a term yields valuable insight into 
the dynamics of more general string models containing non-trivial Ramond-Ramond 
(RR) fields. The interplay between such a cosmological constant and the two-form 
potential had not been considered previously. It was found that the interactions led 
to heteroclinic orbits in the phase space, where the universe underwent a series of 
oscillations between expanding and contracting phases. 

The purpose of the present paper is to extend the work of to both isotropic 
and anisotropic cosmologies containing spatial curvature terms. In particular, we 
consider the class of 'isotropic curvature' universes P,^. These are spatially homoge- 
neous but contain non-trivial curvature and anisotropy. The characteristic feature of 
these models is that the three-dimensional Ricci curvature tensor is isotropic. In other 
words, ^^^Rij is proportional to kSij on the spatial hypersurfaces and these surfaces 
therefore have constant curvature k [§. The class of isotropic curvature universes 
contains the Bianchi type 1 {k = 0) and V (A; < 0) models and a special case of the 
Bianchi type IX (/c > 0) models 0. 

The four-dimensional line element of these cosmologies is given by 



where the one-forms {u^ , u'^ , u^} are given in The spatial hypersurfaces, t = 

constant, are the surfaces of homogeneity. The variable, a, parametrizes the effective 
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spatial volume of the universe and (3 determines the level of anisotropy. We refer to it 
as the shear parameter. In essence, the isotropic curvature models can be regarded as 
the simplest anisotropic generalizations of the flat (A; = 0), open {k < 0) and closed 
{k > 0) FRW universes, respectively. The isotropic models are recovered when (3 = 0. 

The paper is organized as follows. In Section 2, we summarize the form of the effec- 
tive actions we consider. The cosmological field equations are derived and the general 
asymptotic behaviour of the two-form potential and spatial curvature is discussed. 
The global qualitative dynamics for the different classes of models is determined in 
Sections 3 and 4. We summarize and conclude our results in Section 5. The func- 
tional forms of the equilibrium points that arise in the analysis are presented in the 
Appendix. 



2 Cosmological Field Equations 



2.1 String Effective Action 

The NS-NS sector of massless bosonic excitations is common to both the type II and 
heterotic superstring theories. The four-dimensional, string effective action for the 
NS-NS fields can be written as 
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S 



d X\ 



^H^^^H^''^ - 2A 



(2.i: 



where the string coupling, = e*, is determined by the value of the dilaton field, 
the space-time manifold has metric, Qf^i^, and Ricci curvature, R, the antisymmetric 
two-form, Bfj^i,, has a field strength H^^x = d[^B^x], and g = detg^u- The central 
charge deficit of the theory is denoted by the constant, A. The value of this term 
depends on the conformal field theory that is coupled to the string and it can be 
positive or negative. Such a term may also arise from the compactification of higher- 
dimensional form fields or non-perturbative corrections to the self-interaction of the 
dilaton field 
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The three-form H^^x is dual to a one-form in four dimensions and the field equa- 



tion for Bfj_i, is solved by the ansatz |T6| 



0", 



(2.2) 



where e^"'^'^ is the covariantly constant four-form. The scalar a may be interpreted 
as a pseudo-scalar 'axion' field. The Bianchi identity, d[^HyXK] = 0, is then solved by 
reinterpreting this constraint as the field equation for a and this latter equation can 



be derived from the dual effective action |T6| 



= y d^x^/^e' 



i?+(V$)'-ie2*(V(T)'-2A 
2 



(2.3) 
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We establish the dynamics of cosmological models derived from Eq. ( p.3|) in Section 
3. 

In paper II, the action 



S = J d'^Xy/^ |e * 



i?+(V$)'--e'*(Va) 



A 



M 



(2.4) 



was considered, where Am was interpreted as a phenomenological cosmological con- 
stant arising from the interaction potential of a slowly rolling scalar field. We now 
proceed to discuss this action within the context of the massive type IIA supergravity 



theory in ten dimensions |]I0[. This theory represents the low-energy limit of the type 
IIA superstring and has been the subject of renewed interest recently following the 
advances that have been made in our understanding of the non-perturbative features 
of string theory [^,[1^. 



In this theory, the NS-NS two-form potential becomes massive. In the absence 
of such a field, the action is given in the string frame by |T2 



S= / d''xy^\e-' i?io + (V^ 



4 
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(2.5) 

where <^ represents the ten-dimensional dilaton field, is the cosmological constant 
and A = (0, 1, . . . , 9), etc. The antisymmetric field strengths Fab and Fabcd for the 
one-form and three-form potentials represent RR degrees of freedom because they 
do not couple directly to the dilaton field The massless theory is recovered when 
m = 0. 

Maharana and Singh ||lT| have employed the Kaluza-Klein technique |jl8|,|19[ to 
compactify the ten-dimensional theory (|2.5| ) on a six-dimensional torus. We consider 
a truncation of the dimensionally reduced four-dimensional action, where we include 
only the components of the four-form on the four-dimensional external spacetime. We 
therefore neglect the moduli fields arising from the compactification of the form fields 
and the gauge fields originating from the higher-dimensional metric. The effective 
four-dimensional action is then given by 



S 



d Xy 



-9 



i? + (V*)^ - 6 (V7)^ 



-48^ 



67P 



^uXk 



1 

-'I 

2 



(2.6) 



where the four-dimensional dilaton, $, is defined in terms of the ten-dimensional 
dilaton by $ = <^ — 67 and 7 parametrizes the volume of the torus. 

The Bianchi identity for the three-form potential is trivially satisfied in four di- 
mensions and its field equation is solved by the ansatz F'^"^'^ = Qe~^'^e'^'^^'^ , where 
Q is an arbitrary constant and e'^"^'^ is the covariantly constant four-form. Applying 
this duality, together with the conformal transformation 



— 



(2.7) 
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implies that the effective four-dimensional action may be expressed in the Einstein 
frame as 



S 



d xJ—g 



i?-i(v«i>)'-i(vy)'--gv 



1 -2^2<S>-V3y _ 1^2g2$+v^|/ 



(2i 



where y = v I27 represents a rescaled modulus field. 

Without loss of generality, one may perform linear translations on the dilaton and 
modulus fields such that the parameters Q and m become effectively equal. Thus, 
action (p.8|) may be written in the form 



S 



d x\ —g 



Vyf -Kme^^ cosh (VSy) 



(2.9) 



where Am is a positive-definite constant. It follows, therefore, that the four-form and 
ten-dimensional mass parameter together provide an effective potential for the mod- 
ulus field with a global minimum located at ?/ = 0. Consequently, the internal space 
can become stabilized for this specific compactification and a consistent truncation of 
the effective action is therefore given by specifying j/ = in Eq. ( p.9| ). Applying the 
inverse of the conformal transformation (|2.7|) then implies that the effective string 
frame action is given by Eq. (2]^) when the axion field is trivial. Thus, a truncated 
form of Eq. (|2.4| ) is relevant to the type IIA theory. The dynamics of these models is 
considered in Section 4. 

For the purposes of deriving the cosmological field equations from the effective 
actions (|2.3|) and (|2.4|), we combine these two expressions into the single action 



S 



d'^X\/^ \ e 



ie^* [Vaf - 2A 



-A 



M 



(2.10) 



where it is understood that either A or Am should be set to zero. We assume that all 
massless degrees of freedom are constant on the surfaces of homogeneity, t = constant. 
The field equations derived from action (|2.1CI|) for the isotropic curvature metric ( |1.1| ) 
are then given by 



a — aif — -p + K + 


lAMe'^+='" = 





(2.11) 


20 _ 02 _ 3^2 _ g^2 ^ _ 


-3K + 2A = 





(2.12) 









(2.13) 




k + 2ak = 





(2.14) 




p + Gap = 





(2.15) 



together with the generalized Friedmann constraint equation 



3a^ - 02 + 6/52 + -p-3k + 2A + AMe^+=^" = 0, 



(2.16) 
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where 

ip = <^-3a (2.17) 
defines the 'shifted' dilaton field 0,0' 



p = ^262^+*^" (2.18) 
may be interpreted as the effective energy density of the pseudo-scalar axion field 0, 

= 2A;exp(-2a) (2.19) 



represents the spatial curvature term, a and f3 are defined in Eq. (|1.1| ), and a dot 
denotes differentiation with respect to cosmic time, t. 

The moduli fields that may also arise in the string effective action from the com- 
pactification of higher dimensions have not been included in Eq. ( p.l|) and it is as- 
sumed, in particular, that the internal dimensions are fixed. We emphasize, however, 
that in the NS-NS model (Am = 0), one may readily include the dynamical effects of 
these extra dimensions in the case of a toroidal compactification, where the internal 
space has the topology T = x x . . . x , hj reinterpreting the shear parameter, 
P, in the field equations ( |2.11| )-( p.l6D [§|,|T^]- Modulo a trivial rescaling, the moduli 



fields, Prn, that parametrize the radii of the circles, S^, have the same functional form 
in the field equations as the shear term. These moduli may therefore be combined 
with the shear by replacing with an expression of the form = + Pm- 

2.2 Asymptotic Behaviour 

Before concluding this Section, we make some general remarks regarding the asymp- 
totic behaviour of the cosmological models discussed above. 

It can be shown that the action (|2.10D is invariant under a global SL(2, R) trans- 



formation acting on the dilaton and axion fields when the cosmological constants A 
and Am vanish |l^ . The symmetry becomes manifest in the conformally related Ein- 
stein frame (|2.71) and may be employed to generate a non-trivial axion field from a 



solution where such a field is trivial ||2T| . Solutions containing a dynamical axion field 
are known as 'dilaton-axion' cosmologies. The constant axion solutions are referred 
to as 'dilaton-vacuum' solutions and are presented in Eq. (|A.1| ). 



The functional form of the dilaton-axion cosmologies has been derived |^T|. The 
general feature exhibited by these models is that the axion field is dynamically im- 
portant only for a short time interval. The solutions asymptotically approach one 
of the dilaton-vacuum solutions ( [A.l|) in the high and low curvature regimes. The 



axion field also results in a lower bound on the value of the dilaton field and therefore 
the string coupling. However, the SL(2, R) symmetry is broken when either of the 
cosmological constants is present in the action (|2.10|) p2| and analytical solutions are 



not known in this case, even in the isotropic limit {(3 = 0). 
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In addition, the variables K and p in Eqs. ( p.l4|) and ( p.l5| ) may in general be 
combined to define the new variable 

S = P-^. (2.20) 
p + K 

This implies that Eqs. ( p.l4| ) and ( |2.15| ) are equivalent to the evolution equation 



2d(l-S2). (2.21) 



Hence, all of the equilibrium points occur either for d = or for = 1 and this 
implies that either p = (S = — l)ori^' = (S = +l) asymptotically if d 7^ 0. 

In general, Eqs. (p.ll| )- (|2.16| ) define a four-dimensional dynamical system. (Al- 
though there are five ODE's, Eq. (|2.16|) may be employed to globally reduce the 



system by one dimension). In the cases in which all of the equilibrium points lie 
on = 1, the asymptotic properties of the string cosmologies can be determined 
from the dynamics in the three-dimensional sets p = or = 0. The latter three- 
dimensional dynamical system was studied in papers I and II. 

In view of this, we explicitly examine the three-dimensional p = case in what 
follows. The only case in which there exist equilibrium points with d = but 7^ 1 
occurs in the NS-NS case (Am = 0) in which K > Q and A > 0. We shall examine 
the full four-dimensional system in this case, although the three-dimensional subset 
p = still plays a principal role in the asymptotic analysis. 

All solutions represented by the equilibrium points that arise in this work are 
presented in the Appendix. Some of these points also arose in papers I and II, where 
they were labelled differently. In Table |1|, we list all the equilibrium points obtained 
and unify the notation employed in the different works. 



3 Non-Zero Central Charge Deficit (Am = 0) 

In this Section we perform the qualitative analysis for the NS-NS string effective 
action ( p.lO| ) with Am = for an arbitrary central charge deficit. Through Eq. 
( p.l6| ), we eliminate the variable p from the field equations, and make the following 



definitions: 

x.f , ..f, z.f, ..f:, 1.4. (3.) 

The ± signs in the definitions for U and V are to ensure that f/ > and V > 
when necessary. With these definitions, all variables are bounded such that < 
{X^ Z, U,V} <l and Eq. ( pJ6D now reads 



= Y'^ ±U TV - X"^ - Z >{). (3.2) 
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This Paper 


Equation 


Paper I 


Paper II g 


(^1,^2)* 




L±, V, {R, A, S, Si, S2y 


w, w^, {Su,,, Si, S2y 




A3 (A3) 


(cy 




St 


O] 




F 


R, A 


K7 




R, A 










N 










Table 1: The equilibrium points/sets obtained in this paper are listed in the first col- 
umn and their location in the Appendix is given in the second column. The equivalent 
equilibrium points (with different notation) in papers 1 1^ and II are listed in the 
appropriate row. *Note: The equilibrium points in parentheses in the first row are just 
endpoints to the lines L±, VorW. Similarly, equilibrium point C in the second row 
is the endpoint for the line Li for Yi = ¥2 = 1. (See the text for the definition ofYi). 
The points R and A in the first row represent different solutions to those represented 
by the points R and A in the fourth row. Also, in the third row, point F corresponds 
to point Si only. 



The variable ^ is defined in each of the following four cases (the section in which 
they occur is indicated in parentheses) by: 

• A > 

- k > Q: = m + ^'^ (Section 3.1), 

- < 0: ^2 ^ 02 (Section 3.2), 

• A < 

- > 0: ^2 ^ 3i^ + v?^ - 2A (Section 3.3), 

- < 0: ^2 = 02 _ 2A (Section 3.4). 

For example, consider A > with > 0; for this case ¥"^ + 11 = 1 and Eq. ( |3.2| ) 
reads 

^pr^ = l- V -X"^ - Z>Q. (3.3) 

Hence, we use {X, V, Z} as the phase space variables (see Section 3.1 for details). We 
now consider each of these cases in turn and introduce subscripts i = (1,2,3,4) to 
the variables {X, Y, Z, U, V} to distinguish the different subsections. 

For each case, we will set up the four-dimensional dynamical system, followed by a 
discussion of the K = invariant set as examined in Paper I. Note that in Section 3.2 
the K = case is identical to that presented in Section 3.1 and is therefore omitted 
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from that section. Similarly, there will be no i^T = discussion in Section 3.4 since 
it is discussed in Section 3.3. We will then examine the p = invariant set. As 
discussed in Section 2.2, all equilibrium points discussed below which have a = 
also have 7^ 1, and hence nearly all the orbits asymptote towards the equilibrium 
points in one of the invariant sets p = or K = 0. The qualitative behaviour of the 
four-dimensional phase space in each case is then examined. 



3.1 The Case A > 0, A' > 

We define = 0^ + SK and utilize the positive signs for Ui and Vi as defined by Eq. 
From the generalized Friedmann equation we have that 



0<X^ + Zi + Vi< 1, + f/i = 1, 



(3.4) 



and therefore we may eliminate Ui (which is proportional to K), and consider the 
four-dimensional system of ODEs for < {Xf, , Zi,Vi} < 1: 



dXi 

It 

dYi 
If 
dZi 

It 

dVi 

It 



^ (1 - Xl) (2 + Y^) - v^(Zi + V{) + (1 - Xl - 



V3 
l-Y^ 



Xl + + -^X^Yi 



2Zi 
-2Vi 



Y,{i-Xl-Z,) + ^X,{l-Yi 
Y.{xl + Z;)-l=X,(l-Yi 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



The invariant sets Y^ = 1, Xf + Zi + Vi = 1, Zi = and Vi = define the boundary of 
the phase space. The equilibrium sets and their corresponding eigenvalues (denoted 
by A) are 



= ±1, = 1 - Xl = 0; 
(Ai, A2, A3, A4) = [V3 ± Xl] , t2, 0, -2V3 



Xi = 0,^1 = 0,^1 = -(2 + 1^2) 



(A±, A2, A3 
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,2Fi,0 



(3.9) 



(3.10) 



The zero eigenvalues arise because these are all lines of equilibrium points. Here, the 
global sources are the lines Li (for Yi > or > 0) and L~ (for Xi < -^). The 
global sinks are the lines Li (for Yi < 0) and (for Xi > ^). 

This case is different from the other three cases to be considered in this Section, 
since it is the only one with the line of equilibrium points, Li, inside the phase space. 
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This line acts as both sink and source, and corresponds to the exact static solution 
( |A.3| ) which generalizes the static 'linear dilaton-vacuum' solution (|A.4| ) This 
solution was examined in for > and was shown by a perturbation analysis to 
be a late-time attractor. 

The lines correspond to the spatially flat, dilaton-vacuum solutions ( |A.1| ). The 



corresponding stable solutions are in the range > ^ for and < \^ for L 
respectively. 

3.1.1 The Invariant Set ^ = for A > 

This invariant set was studied in paper I, and the dynamics there used the variables 
{Xi, Zi,v = 1 — Xf — Zi — Vi}. It was found that Zi is a monotonically increas- 
ing function (as it is in the full four-dimensional set). The early time behaviour of 
most trajectories is to asymptote towards the linear dilaton-vacuum solution ( |A.4| ), 
represented by the point C, where all degrees of freedom except if are dynamically 
static. To the future, these solutions asymptote towards the line L"*" for Xi > — 4^. 
We note that the point C is the Yi = 1 endpoint of the line Li. Fig. 1 depicts this 
phase space. 

3.1.2 The Invariant Set p = for A > 0, K > 

In the p = case, the system reduces to the three dimensions of {Xi, Yi, Zi} (Vi = 
1 — Xf — Zi). The equilibrium points are the lines with eigenvalues Ai, A2 and 
A3 (from Section 3.1), and the two endpoints of Li with eigenvalues \± and A3 (from 
Section 3.1). These endpoints are specified by the condition = 1 and we denote 
them by Li^\ where the "±" in the superscript reflects the sign of Yi. 

For this invariant set the entire line L"*" acts as a global sink and the entire line 
acts as a global source. Furthermore, we note that for Lf, A_ = 0, and so these 
two points are non-hyperbolic. However, the eigenvectors associated with these zero 
eigenvalues are both 1, 0] and are completely located in the (Xi, Yi) plane. Hence, 
if we choose Zi = and rotate the (Xi, Yi) axes such that 

x^{Y,Tl)- y^(Y,Tl) + ^Xi, 

we see that in the vicinity of the equilibrium point, the trajectories along x for y = 
are also along these eigenvectors. Hence, for y = and small x, it follows that 

dx X 
dT ^7' 

Consequently, for Yi = -|-1, the trajectory along x asymptotes towards the equilib- 
rium point, whereas the trajectory along x for Yi = — 1 asymptotes away from the 
equilibrium point. This implies that the points Lf are saddle points. The phase 
space is depicted in Fig. 2. 



9 



▲ 




Figure 1: Phase diagram of the system ( \3. 3i )-( \37^ ) in the NS-NS (A > 0) sector with 
p ^ and K = 0. The label refers to a line of equilibrium points. In this phase 
space, > is assumed. The labels in all figures correspond to those equilibrium 
points discussed in the text. Throughout, large black dots will represent sources (i.e., 
repellers), large grey- filled dots will represent sinks (i.e., attractors) , and small black 
dots will represent saddles. Grey lines represent typical trajectories found within the 
two-dimensional invariant sets, dashed black lines are those trajectories along the 
intersection of the invariant sets, and solid black lines are typical trajectories within 
the full three-dimensional phase space. 
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The quantity Xi/a/Zi is monotonically decreasing. Such a monotonic function 
excludes the possibihty of periodic or recurrent orbits in this three-dimensional space. 
Therefore, solutions generically asymptote into the past towards L~ and into the 
future towards L"*". In this three-dimensional set, spatial curvature is dynamically 
important only at intermediate times. 



L~ 




Figure 2: Phase diagram of the system 3i )- ^37^ ) in the NS-NS (A> 0) sector with 
p = and K > 0. The labels L"*" and L~ refer to lines of equilibrium points, and the 
labels L^{^^ and l[ ■* represent the equilibrium points which are the endpoints of the 
line Li (for which Xi = +1 and Xi = —1, respectively). In this phase space, > 
is assumed. See also caption to Fig. 1. 

3.1.3 Qualitative Analysis of the Four-Dimensional System 

The qualitative dynamics in the full four-dimensional phase space is as follows. The 
past attractors are the line L~ for Xi < ^ and the line Li for Yi > 0. The future- 
attractor sets are the lines for Xi > ^ and the line Li for Yi < 0. We note 
that Yi/^yVi is monotonically increasing and this implies that there are no periodic 



11 



or recurrent orbits in the full four-dimensional phase space. Therefore, solutions are 
generically asymptotic in the past to either the hne L~ for /i^, < | or to the line Li 
for n > 0. Similarly, solutions are generically asymptotic in the future to either L"*" 
for h^: > ^ or Li for n < 0. 

Since Y\l \fV\ is monotonically increasing, Yi — >■ +1 or Vi asymptotically to 
the future. These limits represent global sinks on the line L+. Conversely, Yi — s> — 1 
or V^i ^ asymptotically to the past, corresponding to the global sources on the line 

. There are also equilibrium points for finite Yi = inside the phase space on 
the line L\. Again, since Y\l is monotonically increasing, the points F^, < are 
global sinks and K,, > are global sources. All of this is consistent with the above 
discussion presented in Section 2.2 regarding the generic asymptotic behaviour. 

We note that the reflections — ^ — Xi and Y\ —Y\ are equivalent to a 
time reversal of the dynamics. Therefore, there are orbits starting on the line L~ 
(for Xx < and ending on the line Li (for Y^ < 0). Similarly, there are orbits 
which begin on the line Li (for Y^ > 0) and end on the line (for Xi > ^). 
Due to the existence of the monotonic function and the continuity of orbits in the 
four- dimensional phase space, solutions cannot start and finish on Li. This is best 
illustrated in the invariant set Zi = 0. In addition, orbits may start on the line L~ 
(for Xl < and end on the line (for Xi > -^). Investigation of the invariant 
set Zi = also indicates which sources and sinks are connected; not all orbits from 
L~ can evolve towards L"*". 

Although the lines lie in both of the invariant sets p = and K = 0, the line Li 
does not. On this line, Xi = 0, and the solutions are therefore static (d = /3 = 0). Eq. 
( ^.21| ) then implies that the axion field and spatial curvature can both be dynamically 
significant at early and late times for the appropriate orbits. 

3.2 The Case A > 0, A' < 

In this case, we choose the negative sign for U2 in Eq. ( |3.1| ), the positive sign for V2, 
and the definition = 0^- The generalized Friedmann constraint (|2.16 ) now implies 
that 

0<Xl + Z2 + U2 + V2<l. (3.11) 

For this system, Y2 = 1, and so the four-dimensional system consists of the variables 

0<{XlZ2,U2,\ 

dX2 
dT 
dZ2 
dT 
dU2 
dT 

dT 



< 1: 






V^(l 


-xl 


-Z2-V2- hj2 


2^2(1 




- Z2) > 0, 


-2U2 


[xl^ 


Z2 + ^X2) , 


-2^2 ( 


xl + 


Z2) < 0. 



X2(l-X2^-Z2j, (3.12) 

(3.13) 
(3.14) 

(3.15) 
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The invariant sets X| + Z2 + U2 + V2 = I, Z2 = 0, V2 = 0, U2 = define the boundary 
of the phase space. The equilibrium sets and their corresponding eigenvalues (denoted 
by A) are 



c 



2 2 4 4 



3' 3' 3' 3 



(Ai, A2, A3, A4) - 

X2 = 0,Z2 = 0,f/2=0,\/2 = l; 

(Ai,A2,A3,A4) = (1,0,2,0), 



Z2 = l-XlU2 = 0,V2 = 0; 



(Ai, A2, A3, A4 



X2 + V3\ , -2, 0, -2^3 



X2 + 



(3.16) 
(3.17) 

. (3.18) 



The point C represents the static 'linear dilaton-vacuum' solution ( [A. 41 ). The 
saddle 5*^ represents the "— " branch of the Milne solution (|A.2|) , where only the 
curvature term and scale factor are dynamic. 



3.2.1 The Invariant Set p = for A > 0, K < 

In the p = case, the system reduces to the three dimensions of {X2, Z2, U2} {V2 = 
1 — X| — Z2 — f/2). The equilibrium points are the same as above with eigenvalues 
Ai, A2 and A3. We note that for this invariant set the entire line acts as a global 
sink, and C acts as a source. Although one of the eigenvalues for the point C is 
zero, it is shown in the following subsubsection that this point is a source in the full 
four-dimensional phase space. The argument is identical in this subsection. 

The variable Z2 is monotonically increasing, and as such we see that the shear 
term is negligible at early times, but becomes dynamically significant at late times. 
There are no periodic or recurrent orbits in this three-dimensional phase space due 
to the existence of this monotonic function. In general, solutions asymptote into the 
past towards the static solution (|A.4|) (point C). Solutions asymptote into the future 
towards L"*". Fig. 3 depicts this three-dimensional phase space. 



3.2.2 Qualitative Analysis of the Four-Dimensional System 

The qualitative behaviour for the invariant set K = is discussed in Section 3.1.1 
(see Fig. 1). In the four- dimensional set, the point C is non-hyperbolic because of 
the two zero eigenvalues. However, it can be shown that this point is a source in 
the four-dimensional set by the following argument. We first note that the variable 
Z2 is monotonically increasing and hence orbits asymptote into the past towards the 
invariant set Z2 = 0. Similarly, V2 is a monotonically decreasing function, and so 
orbits asymptote into the past towards large V2 (i.e., V2 = 1). This implies that they 
asymptote towards the point C. 
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t 



X2 




Figure 3: Phase diagram of the system \3.12j - (^.13[ ) in the NS-NS (K > sector 
with p = and K < 0. Note that represents a line of equilibrium points. See also 
caption to Fig. 1. 
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Since Z2 = asymptotically, let us consider the invariant set Z2 = 0, where Eqs. 
(CT)-(CT) become: 

^ = V3{l-X',-V2-lu2)+X2{l-Xl), (3.19) 
dU2 _ f^2 1 



-2U2 ^Xi + -j=X2 j , (3.20) 
-2^2X2^. (3.21) 



dT V VS 

dT 



It is clear from Eq. (|3.21|) that V2 increases monotonically into the past. Now, this 
three-dimensional phase space is bounded by the surface X| + U2 + V2 = 1, the 
"apex" of which lies at V2 = 1 (and X2 = U2 = 0)- Therefore, all orbits in or on 
this phase space boundary lie below V2 = 1, and therefore asymptote into the past 
towards V2 = 1. To further illustrate that this point is indeed a source, it is helpful 
to consider the invariant set Z2 = U2 = 0, X| + V2 = 1. In the neighbourhood of C, 



Eq. ( p.l9| ) becomes dX2/dT = ^2(1 — X|), indicating that orbits are repelled from 



X2 = 0. Hence, the point C is the past attractor to the full four- dimensional set. 
The future attractor for this set is the line (for X2 > Both C and L"*" lie in 

both of the invariant sets p = and K = 0, which is consistent with the analysis of Eq. 
( p.21| ). We conclude, therefore, that the spatial curvature terms and the axion field 
are dynamically important only at intermediate times, and are negligible at early 
and late times. The dynamical effect of the shear becomes increasingly important 
because the variable Z2 increases monotonically. On the other hand, the variable 
V2 decreases monotonically and the dynamical effect of the central charge deficit. A, 
becomes increasingly negligible. In addition, the existence of monotone functions in 
the four-dimensional phase space prohibits closed orbits and serves as proof of the 
evolution that is described above. 

3.3 The Case A < 0, if > 

We choose the positive sign for f/3 and the negative sign for V3 in Eq. ( |3.1| ) to ensure 
that these variables are positive definite. We also define = 0^ + 3K — 2A for this 
case. The generalized Friedmann constraint equation can then be rewritten as 

0<Xl + Zs< 1, Yi + f/3 + 1/3 = 1, (3.22) 

We again eliminate f/3 (which is proportional to K), and consider the four- dimensional 
system of ODEs for < If, Z3, ^3} < 1: 



dX. 



1-Xi- Z,)[V3 + XsYs) -^{l-Yi- V3){l - Xl - Z,) ,(3.23) 



dT V ^ A ^ V3 

dYs ^ J_ 
dT 73 



"^^^ ),X,Y, (1 - Yi - V,) + (1 - ^3^) (X| + Z3) , (3.24) 
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dT 
dT 



2Z, 



2V. 



J_X,{l~Y,'-Vs)+Ys{l 



.V3 
1 



X3 - Z3 



(3.25) 
(3.26) 



The invariant sets Y^ + Vs = 1, X| + Z3 = 1, Z3 = define the boundary of the phase 
space. The equilibrium sets and their corresponding eigenvalues (denoted by A) are 

L±: Y. = ±hZ. = 1-XlV. 



^3 = ±1,^3 = 1 
(Ai, A2, A3, A4) = 



^3' "3 

2 

73 



0; 



VSiXgl ,^2, 0,^2^3 



±X3 



(3.27) 



where again the zero eigenvalues arise because these are all lines of equilibrium points. 
Here, the global sink is the line for X3 > ^ (saddle otherwise), and the global 
source is the line L~ for X3 < ^ (saddle otherwise). Stable solutions on the line L"*" 
correspond to the range /i* > The stable solutions on L~ arise when < |. 



3.3.1 The Invariant Set fsT = for A < 

This invariant set was studied in paper I, the dynamics of which are as follows. The 
variables X3 and Y^ are monotonically increasing functions, corresponding to a and 
(respectively). The former implies that these trajectories represent cosmologies that 
are initially contracting and then reexpand. In paper I, the third variable used was 
K = 1 — X| — Z3. This is proportional to p and is only dynamically significant at 
intermediate times. It asymptoted to zero into the past and future, indicating that 
the axion field is negligible at early and late times. All the equilibrium points in this 
invariant set are represented by the dilaton- vacuum solutions ( [A.l|) (the lines L^)- 
In general, orbits asymptote into the past towards the line L~ (for X3 < -^), and 
asymptote to the future towards the line (for X3 > ^). Fig. 4 depicts this 
three-dimensional phase space, using the variables {X3,y3,K}. 

3.3.2 The Invariant Set p = for A < 0, ^ > 

Since d 7^ at the equilibrium points, we examine the p = case, where the system 
reduces to the three dimensions {X3, Y^, V3} (Z3 = 1 — X|). The equilibrium points 
are the same as above with eigenvalues Ai, A2, A3. We note that the entire line L"*" acts 
as a global sink and that the entire line acts as a global source in this invariant 
set. The function Ys/y/V^ is monotonically increasing, eliminating the possibility of 
periodic orbits. Fig. 5 depicts this three-dimensional phase space. 



3.3.3 Qualitative Analysis of the Four-Dimensional System 

The qualitative dynamics in the full four-dimensional phase space is as follows. The 
global repellers and attractors are the lines L~ and L+, respectively. Orbits generically 
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K 




Figure 4: Phase portrait of the system (\3.2d^ )- ^^2BJ in the NS-NS fA < Oj sector for 
K = and p 7^ 0. Note that the labels and L~ refer to lines of equilibrium points. 
See also caption to Fig. 1. 



17 




18 



asymptote into the past towards the hne L~ (for X3 < and into the future 

towards L"*" (for X3 > ^). In the former case the stable solutions are given by the 
"+' branch of Eq. (|A.1| ) for /;,* < |. The stable sinks are given by the "— " branch 
with h^: > Again, we see that the curvature term, axion field and the central 
charge deficit are dynamically important only at intermediate times. The existence 
of the monotonically increasing function Is/v^ excludes the possibility of periodic 
orbits and serves to verify the above description of the evolution of the solutions in 
the four-dimensional set. 



3.4 The Case A < 0, < 



For this case, the appropriate definition for the variable ^ is 



2A and we 



choose the negative signs for both U4 and V4 in Eq. (|3.1| ). The generalized Friedmann 
constraint equation is written as 



<Xf + Z. + Ua< 1, 



^4^^4-1, (3.28) 

Treating V4 as the extraneous variable results in the four-dimensional system consist- 
ing of the variables < Y^, Z4, U4} < 1: 



dT 
dY4 

'dT 

dZi 

It 

dUi 

It 



1 



ZA ( + X4F4 



1 - Yi) (xl + Z4) > 0, 

2.^4 y4 ^1 — xl — Z/^ 



4 5 



-2Ua 



I4 + Zij + —j=Xi 



(3.29) 

(3.30) 
(3.31) 
(3.32) 



The invariant sets X| + Z4 + f/4 = 1, Z4 = 0, Y^ = 1,^/4 = define the boundary 
of the phase space. The equilibrium sets and their corresponding eigenvalues (denoted 
by A) are 

1 „ _ 2 

3' 



X, 



±1,^4 = 0,^/4 



(Ai, A2, A3, A4) — 

F4 = ±1,^4 = 1 

(Ai, A2, A3, A4) = 



2 2 4 4 

±-,T-,±-,±- 
3' 3' 3' 3 

XlU, = 0; 



V3±xA,t2,0,t2V3 



73 



(3.33) 



(3.34) 



The global source for this system is the line L~ (for X4 < /i* < |) and the global 
sink is the line (for X4 > h^.. > ^). The saddle points, 5*^, are represented 
by the Milne models ( |A.2| ), where corresponds to the "— " solution and S" to the 
"+" solution. 
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3.4.1 The Invariant Set p = for A < 0, K < 



— Z4) 



This system reduces to the three dimensions of {X4, F4, Z4} {U4 = 1 
The equilibrium points are the same as above with eigenvalues Ai, A2, A3. The 
entire lines and L~ now act as a global sink and source, respectively. Recurrent 
orbits are forbidden by the existence of the monotonically increasing variable Y4. 
Hence, solutions generically asymptote into the past (future) towards the "+" ("— ") 
dilaton-vacuum solutions (|A.1|) and the curvature term and central charge deficit 
are dynamically significant only at intermediate times. Fig. 6 depicts this three- 
dimensional phase space. 




Figure 6: Phase portrait of the system P^^ - P^^^ in the NS-NS (A < 0) sector 
with p = and K < 0. Note that and L~ represent lines of equilibrium points. 
See also caption to Fig. 1. 



3.4.2 Qualitative Analysis of the Four-Dimensional System 

The dynamical behaviour in the invariant set i^' = is identical to that described 
in Section 3.3.1 (see Fig. 4). The qualitative dynamics in the full four-dimensional 
phase space is as follows. Since Y4 is monotonically increasing, the orbits asymptote 
into the past towards Y4 = —1 and into the future towards Y4 = +1. As in the 
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above examples, the existence of such a monotonic function excludes the possibility 
of periodic orbits in the four-dimensional phase space. Most orbits asymptote into 
the past towards the line L~ (for X4 < and into the future towards the line L"*" 
(for X4 > ^). The range of values for h^, corresponding to stable solutions is given 
by > ^ for those on and < | for those on L~ . The variable I4 increases 
monotonically along orbits and this implies that (p is dynamically significant at early 
and late times, but dynamically insignificant at intermediate times. Conversely, we 
see that the curvature term, central charge deficit and axion field are dynamically 
important only at intermediate times, and are negligible at early and late times. 

This concludes the qualitative analysis of the isotropic curvature string cosmolo- 
gies with NS-NS fields. In the next Section, we consider the effect on the dynamics 
of introducing a non-trivial Am. 



4 Non-Zero Cosmological Constant Am (A = 0) 

We begin this Section by defining new variables 

4 = 6^(^+3-) A, AT = 6/3", ^ = v\ h = a', K = ke-^^+'''\ (4.1) 
at dT 

where a prime denotes differentiation with respect to the new time coordinate, T. 
Eqs. ( |2.11|) -( p.l6|) may then be written as the set of ODEs: 



h' 




1 , 1 


2 2 ^ 


2 + ^pe 


^' 






2 2^2 


2 A'^ 


N' 


= {'ilj-3h)N 




K' 


= -{tlj + 5h)K 




P' 


= —6hp 




3h^ 


- ^^ + N-3K + Am 



(4.2) 

(4.3) 

(4.4) 
(4.5) 
(4.6) 

(4.7) 



The variable p may be eliminated due to the constraint equation ( p.l6|) . It also 
proves convenient to further define a set of variables 

V3h ^ N ±3K , ±Am d d 

P^^, X=^, - = ^, C = ^, A = ^, ^=e^, (4.8) 

where the ± signs ensure that C > ^-^^d A > 0. The variable is defined in each of 
the following four subsections by 

• Am > 
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K > 0: ^2 ^ 3ir + (Section 4.1), 
K <0: ^^ = iIj^ (Section 4.2), 



• Am < 

- K > 0: ^'^ = 3K + - Am (Section 4.3), 

- K <0: = iIj^ -Am (Section 4.4). 

With these definitions, all variables are bounded, < {/i^, x^, J^, C; A} < 1, and Eq. 
( [4.7|) now reads 



1 



-2„-(v+3a) 



x'±CTA-/i'-z/>0. 



(4.9) 



Our overall approach in this Section is identical to that of Section 3 and we refer the 
reader to the discussion immediately after Eq. (|3.2|) for the general outline adopted. 



4.1 The Case Am > 0, A' > 

For K > 0, Eq. ( [4. 7] ) is written in the new variables as 

< /i? + z/i + Ai < 1, Ci + X? = 1, 



(4.10) 



where the "+" sign is chosen for both A and ( in Eq. ( [4.8|) . For this case, the variable 
(i = 1 — Xi will be considered extraneous and the system (^.2|)-(^3|) then reduces to 
the four- dimensional system: 



dfii 

dxi 
dr 
dvi 
dr 
dXi 



^1 -^Ai) (v^ - ^ (i (i 



+ + o Ai + —f^l^iXi 



^1 



Ai 



2 

■ 1 
73 



2 "^ ■ VS 
(l - xl) + 2X1 {l-f^l- 

/ii (5 - 2x1) + Xl (l - 2fil - 2z/i - Ai 



^1 - 2A1 



y3Ai,(4.11) 
(4.12) 
(4.13) 
(4.14) 



The invariant sets jj^j + ui + Xi = 1 {p = 0), xl = I {K = 0), i/^ = (A^ = 0) and 
Ai = (Am = 0) define the boundaries to the phase space. The equilibrium points 
and their respective eigenvalues (denoted by A) are given by 



Xl = ±1,/^? + z/i = 1, Ai = 0; 



(Ai, A2, A3, A4 



T2V3 



(4.15) 
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± . 



(Ai, A2, A3, A4 



Xi 




V23T 




(4.16) 



From the eigenvalues, we deduce that is a late-time attractor for /i^ < |, and that 
L~ is an early-time repeller for /i^ < |. In both cases, Xi = 1 and therefore K = 
(Ci = 0). The "+" solution of (^) is a sink for K < ^ and the "-" solution of ( |0| ) 
is a source for < \. The points are saddle points on the boundary of the phase 



space and correspond to the exact solution ( [A .51 ). 



4.1.1 The Invariant Set fsT = for Am > 

This invariant set was studied in paper II by employing dynamical variables equiv- 
alent to {/ii, Ai, i^i}, and we now summarize the important features of the model. 
For spatially isotropic solutions confined to the invariant set i^i = 0, most trajec- 
tories evolve from the equilibrium point Si (labelled "F" in paper II) located at 
(yUi,z/i,Ai) = ("^)0,^) (corresponding to the "— " solution ( |A.5|) ) and are future 
asymptotic to a heteroclinic orbit. There are two saddle equilibrium points, 5*1 and 
5*2, which are the endpoints of the line L"*"; 5*1 is given by the "— " branch of Eq. ( |A.1| ) 
with ^* = ^ and 5*2 corresponds to = There are also the single boundary 
orbits in the invariant sets Ai = 0, corresponding to Am = and Ai + /i^ = 1 (constant 
axion field). An orbit spends the majority of its time in the neighbourhoods of 5*1 
and 5*2 and shadows the respective boundary orbits as it rapidly moves between the 
two saddles. Progressively more time is spent near the saddles for each completed 
cycle, and the dynamics is therefore not periodic. 

An anisotropic contribution (z/i 7^ 0) does not introduce new sources into the 
system and the point Si is still the only source. Fig. 7 depicts the full three- 
dimensional space. Eq. ( |4.13| ) implies that Ui is a monotonically increasing function 
and consequently the orbits are repelled from Si and spiral out monotonically. The 
general behaviour for most trajectories in this phase space is to evolve away from the 
equilibrium point Si, spiral about the line /ii = ui = y(Ai — ||) and eventually 

asymptote towards the line for < |. These trajectories were discussed in detail 
in paper II. 



4.1.2 The Invariant Set p = for Am > 0, fsT > 

For the invariant set p = 0, the system ( |4.11|) -( p:.14|) reduces to the three dimensions 
{/ii, X15 1^1} (Ai = 1 — fii — ui). In this invariant set, /i is a monotonically decreasing 
function, and so the possibility of periodic orbits is excluded. The only equilibrium 
points are the lines with the first three eigenvalues in Eq. (|4.15|) , and therefore 
the early and late time attractors are the lines L~ (for pi > > ^) and (for 
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4 




V. 



1 



Figure 7: Phase portrait of the system ^.ll\ )-( ^.ldi ) for Am > with p 7^ and 
K = 0. Note that the label refers to a line of equilibrium points. In this phase 
space, > zs assumed. See also caption to Fig. 1. 
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/ii < h^. < |). The curvature and cosmo logical constant are only dynamically 
significant at intermediate times and the phase space is depicted in Fig. 8. 




Figure 8: Phase portrait of the system ^.11\ }- ^.13^ } for Am > with p = field and 
K > 0. Note that and L~ represent lines of equilibrium points. See also capt 
to Fiq. 1. 
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4.1.3 Qualitative Analysis of the Four-Dimensional System 

The qualitative dynamics in the full four-dimensional phase space is as follows. The 

1 1 _ 1 

function Xi'^i ~ lA^ ~ ~ -^i) ^ monotonically increases and so there can be no 
periodic orbits. The only past-attractors belong to the line for /i^ < | {hi < |) 
and the only future attractors belong to the hne for /i^ < | [hi < |). Both hues 

lie in both of the invariant sets p = and K = which is consistent with the 
analysis of Eq. ( p.21|) . This implies that the spatial curvature, cosmological constant 
and axion field are only dynamically significant at intermediate times. 

For orbits in the four- dimensional phase space, the complex eigenvalues of the sad- 
dle point Sf suggest that the heteroclinic sequences may exist in the four-dimensional 
set. Indeed, those orbits which asymptote to the K = {] invariant set do generically 
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end in a heteroclinic sequence, interpolated between two equilibrium points represent- 
ing two dilaton-vacuum solutions, as discussed in Section 4.1.1. However, for those 
orbits which asymptote towards the p = invariant set, there are no heteroclinic 
sequences (as is evident from Fig. 8). 



4.2 The Case Am > 0, A' < 

For K < 0, Eq. ( [4.7|) is written in the new variables as 

0</x^ + z/2 + C2 + A2<l, X2 = l, (4.17) 

where the "+" sign for A and the "— " sign for ( have been chosen in Eq. ( |4.8| ). For 
this case, we explicitly choose X2 = +1, as X2 = — 1 corresponds to a time reversal of 
Eqs. ([4.2D-(^^. The system then reduces to the four-dimensional system: 



dfX2 

dr 

dr 
d\2 

dr 



1-/^2 



1 



A2) fV3 



1^2 



2U2{1 - ^2-1^2- -A2 



— = -2C2 l^fil + 1^2+ 2^2 + ^/i2 j 

A2 (1 - 2/i2 - 2z/2 - A2 + VSlJ,2 



VS{\2 + -C2 



(4.18) 
(4.19) 
(4.20) 
(4.21) 



The invariant sets /i^ + + C2 + A2 = 1 (p = 0), C2 = (TT = 0), z/i = (A^ = 0) and 
Ai = (Am = 0) define the boundaries to the phase space. The equilibrium points 
and their respective eigenvalues (denoted by A) are given by 



5+ : 



/^2 



N : 



(Ai, A2, A3, A4) 



1^2 



-,Z^2 



A2 = 0,1^2 = 0; 



3(2,4,-2,4) 



(Ai, A2, A3, A4) 

P2 + ^^2 = 1, C2 

(Ai, A2, A3, A4) 



-(1+2^2,1 

5 

^0,A2 = 0; 



^72, 4, 2) 
P2 + , V3 



(4.22) 



(4.23) 



/^2 - 



,-2v^ 



fJ'2 + 



(4.24) 

From the eigenvalues, it is clear that is a late-time attractor for P2 < | < |)- 
The point inside the phase space is the early-time attractor for the system, and 
represents the curvature-driven solution ( [A.6| ). The saddle point 5*+ corresponds to 
the "-" branch of the Milne solution 
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4.2.1 The Invariant Set p = for Am > 0, K < 



For this invariant set, the four-dimensional system ([4.18|) -( p:.21| ) reduces to a three- 



dimensional system involving the coordinates {/i2, 1^2,^2} (C2 = 1 — pi — 1^2 — A2). 
The equilibrium points are the same as for the full four- dimensional set, but the 
eigenvalues are now (Ai,A2,A3). The variable z/ is a monotonically increasing func- 
tion, the existence of which eliminates the possibility of recurrent orbits. Thus, the 
generic behaviour of this model is for solutions to asymptote into the past towards the 
curvature-dominated solution ( [A .61 ), represented by the point N, and to the future 
towards the line L+ for p2 < {K < |)- Fig- 9 depicts this phase space. 




Figure 9: Phase portrait of the system ^.1^ )-( ^.19{ ) for Am > with p = and 
K < 0. Note that represents a line of equilibrium points. In this phase space, 
if > is assumed. See also caption to Fig. 1. 

4.2.2 Qualitative Analysis of the Four-Dimensional System 

The invariant set K = was described in Section 4.1.1 (see Fig. 7). In the full 
four- dimensional set, the point is the early-time attractor, and L"*" is the late-time 
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attractor for /i^ < | {hi < |). The point hes in the invariant set p = and hes 
in both of the invariant sets p = and K = 0. This is consistent with the analysis of 
Eq. (|2.21|) . We see that z/2 is a monotonically increasing function, and so there are no 



recurrent or periodic orbits. Furthermore, since z/2 increases monotonically, the shear 
in the model is initially dynamically trivial, but becomes significant asymptotically 
into the future. The axion field and cosmological constant are only dynamically 
important at intermediate times and do not play a role in the early- and late-time 
behaviour of the cosmologies. The curvature term is dynamically significant at early 
and intermediate times, but becomes dynamically trivial at late times. 

Orbits which asymptote into the future towards the K = invariant set generically 
end in a heteroclinic sequence, as described in Section 4.1.1 and depicted in Fig. 
7. However, such a sequence does not occur for orbits which asymptote into the 
future towards the p = invariant set. Indeed, by examining the eigenvalues of the 
equilibrium points of the four-dimensional system, there do not seem to be heteroclinic 
sequences outside of the K = invariant set. 



4.3 The Case Am < 0, A' > 

For completeness we now consider the cases where Am < 0. When > 0, Eq. ( |4.7| ) 
is written in terms of the new variables as 



< Ps + ^^3 < 1, X3 + C3 + A3 



(4.25) 



where the "— " sign for A and the sign for ( have been chosen in Eq. 
The variable A3 is chosen as the extraneous variable and the system ( [4.2|) -( ^l5D then 
reduces to the four-dimensional system: 



dr 
dr 



{1- fil- Us) {VS + P3X3) + ^ (1 - P3 



I-Xl- ICs) , (4.26) 



2 



/i3X3 (1 - X3 - ^(3)] - ^ (1 - xl) (1 - - 2us) + ^Cs 



dP3 
dr 

dr 



2X3(l-/i3-'^3 



(i-xl 



^ = -C3 



2X3 {f^l + '^3) + ^/i3 (5 - 3x3 - 5Cc 



(4.27) 
(4.28) 

(4.29) 



The invariant sets p| + 2/3 = 1, xl + C3 = 1? ^^3 = and C3 = define the boundary to 
the phase space. The equilibrium sets and their corresponding eigenvalues (denoted 
by A) are 



X3 = ±l,/i3 + ^^3 = 1,(3 = 0; 
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R : 



A : 



(Ai, A2, A3, X4 



Vs. 



X3 



V3 

(Ai, A2, A3, A4 
1 



,fl3 = -l,Z/3 = 0,C3 



0; 



-2V3 



^(1,2,6,10) 



X3 = ' ^3 = 1, = 0, Cs = 0; 



(Ai, A2, A3, A4 



1 



1,2,6,10). 



(4.30) 



(4.31) 



(4.32) 



Here, there are two early-time attractors. The first is the point R, representing the 
"+" branch of the constant axion, spatially isotropic solution ([A.7|) , where (p < 0. 
The second is the line L~ for < | (hi < |). Likewise, there are two late-time 
attractors: these are the point A, representing the "— " solution of Eq. ( |A.7| ), where 
> 0, and the line L+ for /ii < | {hi < |). 



4.3.1 The Invariant Set K = for Am < 

In paper II, this invariant set was examined using variables which are the same as 
{/i3, Xs, i^s}- There it was shown that /is is monotonically increasing, so that most 
trajectories in this phase space represent bouncing cosmologies which are initially 
contracting. Most trajectories asymptote into the past towards either L~ or to R 
(see above). To the future, most orbits asymptote towards either or A. Fig. 10 
depicts this three-dimensional phase space. 



4.3.2 The Invariant Set p = for Am < 0, fsT > 

For this invariant set, the four-dimensional system ([4.26|) -( p:.29|) reduces to a three- 
dimensional system involving the coordinates {fJ^SyXsXs} (^3 = 1 ~ A*!)- The equi- 
librium points are the same as the full four-dimensional set, but with eigenvalues 
(Ai, A2, A3), and so the line L~ is a source for ^^3 > ^ (^* > ^) and the line L"*" is a 
sink for A^s < {h* < |). The function Xs/y^ is monotonically increasing, and so 
there are no recurring or periodic orbits. Hence, solutions generically asymptote into 
the past towards L~ or R and into the future towards L"*" or A. 



4.3.3 Qualitative Analysis of the Four-Dimensional System 

The function (^X3 + ^/^s) monotonically increasing, and so there are no re- 

curring or periodic orbits in the full four- dimensional set. The early-time attractors 
are the line for < | {h^ < |) and the point R. The late-time attractors are 
the line for /i| < | (hi < |) and the point A. All of these attractors lie in 
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^^3 




Figure 10: Phase portrait of the system ( 14-26^ )- ^. 2!^ ) for Am < with p ^ and 
K = 0. Note that the labels and L~ refer to lines of equilibrium points. See also 
caption to Fig. 1. 
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Figure 11: Phase portrait of the system ( \4.2di )- ^.2[^ for Am < with p = and 
K > 0. Note that and L~ represent lines of equilibrium points. See also caption 
to Fig. 1. 
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both of the invariant sets p = and K = and the axion field and curvature term 
are dynamically significant only at intermediate times. For early and late times, the 
cosmological constant is dynamically important only when the shear is dynamically 
trivial (e.g., the points R and A), and vice- versa (e.g., the lines L^)- 



4.4 The Case Am < 0, A' < 

For this case, Eq. 



M ^ 

is written in terms of the new variables as 



< pa + z/4 + C4 < 1, 



xl + ^4 



(4.33) 



where the "— " sign for both A and ( has been chosen in Eq. ([4.8| ). The variable A4 
is chosen as the extraneous variable and the system ([4.2|)-(P3|) then reduces to the 
four-dimensional system: 



(ip4 
dr 

dr 
(iz/4 

dr 
dU 
dr 



{l- z/4 
1 

~2 



V3 + 



P4X4 



+ ^ (1 - (1 - - ^C4, (4.34) 



1 - xt) fl - 2p4 - 2z/4 + V3p4X4 



z/4 

-C4 



2x4 (1 - P4 - z/4) - V3p4 (l - xl 



2X4 (P4 + ^i) + (5 - 3X4 



(4.35) 
(4.36) 
(4.37) 



The invariant sets P4 + z/4 + ^4 = 1, xl = 1, ^^4 = and C4 = define the boundary to 
the phase space. The equilibrium sets and their corresponding eigenvalues (denoted 
by A) are 

: X4 = ±1,/U4 + ^^4 = 1,(4 = 0; 

2 



(Ai, A2, A3, A4 



0,T 



-\/3 ± U4 , -\/3 



P4 =F 



%/3 



-2V3 



P4 ± 



1 
73 



R : 



X4 = = -l,z^4 = 0,(4 = 0; 



^(1,2,6,10) 



A 



(Ai, A2, A3, A4) 

" "71' = ^' ^4 = 0, C4 = 0; 



(Ai, A2, A3, A4 



1 

73 



(1,2,6,10), 



^ 2 

X4 = ±i,P4 = "^'^^4 = 0,(4 = -; 
(Ai,A2,A3,A4) = ^(Tl,±l,±2,0). 



(4.38) 



(4.39) 



(4.40) 



(4.41) 
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As in the previous case, there are two early-time attractors: the point R and the hne 
L~ for /i| < I (hi < |). Likewise, there are two late-time attractors: the point A 
and the line for /u| < | (hi < |). 



4.4.1 The Invariant Set p = for Am < 0, K < 



For this invariant set, the four-dimensional system ([4.34|) -( p:.37|) reduces to a three- 



dimensional system involving the coordinates {fi^, Xi^ ^i} id = ^ ~ f^l ~ ^i)- The 
equilibrium points are the same as the full four- dimensional set, and the eigenvalues 
are (Ai, A2, A3). The difference is that the line L"*" is a sink for /i4 < ^ and the line 
L~ is a source for /i4 > The function ^4^/ ^/Ul is monotonically increasing, and 
so periodic orbits cannot occur. Hence, solutions generically asymptote into the past 



towards the "+" branch of Eq. ( [OD for K> ^,01 the solution of Eq. 



Into the future, solutions asymptote towards either the "— " branch of Eq. ( |A.1| ) for 
< |, or to the "— " solution of Eq. (|A.7| ). Fig. 12 depicts this phase space. 




Figure 12: Phase portrait of the system (^T^j-(^75^ for Am < with p = and 
K < 0. Note that and L~ represent line of equilibrium points. See also caption to 
Fig. 1. 
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4.4.2 Qualitative Analysis of the Four-Dimensional System 

The invariant set = is discussed in Section 4.3.1 (see Fig. 10). The function 
lii/ s/iM is monotonically increasing, and so there are no recurring or periodic orbits. 
The asymptotic behaviour is identical to that described in Section 4.3.3. 



5 Discussion 

In this paper, we have presented a complete qualitative analysis for the isotropic 
curvature string cosmologies derived from the effective action ( ^.lUD in the two cases 
where either A = or Am = (the cases where both terms are non-zero is discussed in 
P^). This was made possible by compactifying the phase space in terms of suitably 
defined variables. When Am = 0, Eq. ( |2.10|) represents the action for the NS-NS 
fields that arise in both the type II and heterotic string theories when an arbitrary 
central charge deficit is present. We identified the cosmological constant Am with 
terms that arise in the RR sector of the massive type IIA supergravity theory when 
appropriate conditions apply. 

The subset /5 = corresponds to the class of spatially isotropic FRW universes 
with arbitrary spatial curvature. In the positively curved case, we have extended the 
work of Easter et al. 0, who performed a perturbation analysis on the static, closed 
FRW model to show that it was a late-time attractor. More generally, the models we 
have considered represent Bianchi type I, V and IX universes. 

For each case, we have established the existence of monotonic functions which 
precludes the existence of recurrent or periodic orbits. Consequently, the early- and 
late-time behaviour of these models can be determined by analysing the nature of 
the equilibrium points/lines of the system. In all cases, the spatially fiat dilaton- 
vacuum solutions L^, given by Eq. ( |A.1| ), act as either early- or late-time attractors 
and, in many cases, act as both. Because these solutions lie in both the p = and 
K = Q invariant sets and contain neither a central charge deficit nor a non-zero 
Am contribution, we may conclude that the shear and dilaton fields are dynamically 
dominant asymptotically. Furthermore, with the exception of the A > 0, iiT > case, 
all early-time and late-time attracting sets lie in either the p = invariant set or the 
K = Q invariant set, and a majority of these sets lie in both. 

Thus, we see a generic feature in which the curvature terms and the axion field 
are dynamically significant at intermediate times and are asymptotically negligible 
at early and late times. The exception to this generic behaviour is the A > 0, i^' > 
case, where the generalized linear dilaton-vacuum solution ( [A.3| ), in which neither 



p = nor K = acts both as a repeller (for < 0) and as an attractor (for > 0). 
In these solutions, the variables p and K are proportional to the central charge deficit, 
A. Note that asymptotically d = (and /3 = 0) and hence these models are static. 

When A < 0, the central charge deficit is dynamically significant only at interme- 
diate times, and is asymptotically negligible at early and late times. In fact, the only 
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repelling and attracting sets in this instance are the dilaton-vacuum solutions. When 
A > 0, the central charge deficit can be dynamically significant at both early and 
late times, and the corresponding solution is the generalized linear dilaton-vacuum 
solution (|A.3| ), represented by the line Li. When > 0, these solutions can be 
repelling (0 < 0) and attracting (0 > 0). When K < 0, the endpoint of this line, C 
(representing Eq. ( [A .41) ) is a repeller. 

When Am > 0, the cosmological constant may play a significant role in the early 
and late time dynamics. For instance, although in the four-dimensional sets there are 
no repelling or attracting equilibrium points in which Am is dynamically significant, 
we have found that the orbits which are attracted to the K = invariant set end in a 
heteroclinic sequence which interpolates between two dilaton-vacuum solutions (see 
Section 4.1.1, Fig. 7 and paper II). During this interpolation, the orbits repeatedly 
spend time in a region of phase space in which Am is dynamically significant (the 
region Ai > in Fig. 7), although most time is spent near the dilaton-vacuum saddle 
points where Am is dynamically negligible. When Am < 0, the cosmological constant 
can be dynamically significant at both early and late times, since solutions typically 
asymptote to the solution (|A.7| ), where the shear and axion field are static. For the 
repelling and attracting sets of the Am < cases, the shear term is only dynamically 
significant when the cosmological constant is not, and vice versa. 

The parameter f3 measures the degree of anisotropy in the models. If we define 
isotropization by the condition jS —>■ [^, then we note that in general f3 ^ 
at the equilibrium points on the line L^. Therefore, solutions asymptoting towards 
the sinks on these lines do not isotropize to the future. At all other equilibrium 
points, however, /3 = and the corresponding appropriate string cosmologies therefore 
'isotropize'. This is an important result and a similar situation occurs in the time- 
reverse models. On the other hand, the question of isotropization of string cosmologies 
in a more general context remains an open question. Note that the shear in the models 
that we have discussed is essentially of Bianchi type I. In general relativity with a 
perfect fiuid, it is known that Bianchi type I models isotropize whereas in general 
spatially-homogeneous models do not isotropize [25, 26]). 

For the equilibrium points corresponding to the sinks on the line L"*", jS"^ oc 0^ 
{p = K = 0), and so the energy densities of the modulus and dilaton fields are 
proportional to one another. Hence, the corresponding dilaton-vacuum solutions are 
"matter scaling" string cosmology solutions, which act as local attractors, similar 



to the matter scaling solutions in general relativistic scalar field cosmologies [27 



Finally, for every equilibrium point within these phase spaces, the scale factor of the 
corresponding exact solution is a power-law function of cosmic time, and therefore 



all of the corresponding exact solutions are self-similar |2^ . 

The cases in which ip < are related to a time reversal of the models discussed in 
the text. We have not explicitly considered these models here; however, the conclu- 
sions concerning isotropization are similar although the details of these models may 
be different (e.g., the role of sources and sinks can be interchanged). 
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In conclusion, therefore, we have estabhshed the quahtative properties of all the 
isotropic curvature string models discussed in the text by finding appropriate mono- 
tone functions. Typically, the curvature term is dynamically significant only at in- 
termediate times and is asymptotically negligible. There are only two exceptions to 
this. The first corresponds to the case {A> 0,K > 0}, where the generalized linear 
dilaton-vacuum attractors and repellers have a non-negligible curvature. The second 
case is the {Am > 0, < 0} model in which the repeller N represents the curvature- 
driven solution ( |A.6D . Finally, we note that when Am > there exist heteroclinic 
sequences in the invariant set K = 0. This implies that the qualitative behaviour as- 
sociated with heteroclinic sequences is only valid for solutions which approach K = 0. 
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A Equilibrium Points 



In this Appendix, we present the analytical solutions to Eqs. ( p.ll|) -( p.l6| ) that 
represent all of the equilibrium points that arise in this paper. 

The 'dilaton-vacuum' solutions correspond to solutions where the axion field is 
constant and cosmological constants vanish (A = Am = o" = 0). In the spatially flat 
case {K = 0), they are power laws: 



a* \t 



e = e \t\ 



g/3 ^ |^|±V(l-3h2)/6 ^ 

k = 0, (A.l) 

where {a*, a*, /;,*} are constants, a = e" is the averaged scale factor of the 
universe, the ± sign corresponds to the sign of t and e = ±1. These solutions have 
a curvature singularity at t = 0. Note that the shifted dilaton field ( p. IT)) satisfies 
if > for t < and < for t > 0. The "— " branch of Eq. (|A.1|) corresponds to 
the line L"*" throughout this paper and the branch corresponds to the line L~ . 
Another solution which appears when both A = Am = and cr = is the Milne 
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form of flat space: 



a = a* (±t) , 

$ = 

13 = 13., 

(J = a^, 

k = -al (A.2) 

where {a*, cr*} are constants. The "±" sign corresponds to the sign of t. These 

solutions are labelled S"^ throughout the paper (note 5"+ corresponds to t < and 
S~ corresponds to t > 0) and arise as saddle points. 

There is a line of equilibrium points that arises when the central charge deficit, 
A, is included in the action ( |2.1CI| ). This class of solutions has the form 



a = 




$ = 




/9 = 




a = 




k = 




2 + n2^ 




2 f 1 - n2 



■ exp — — n 




Aa^, (A.3) 

where n G [—1, 1]. The solution is static, but has non-trivial spatial curvature, dilaton 
and axion fields. These solutions are represented by the line Li in Section 3.1. In 
that Section they represent past attractors for n G (0, 1] and future attractors for 
n G [—1,0). These solutions were found in 0] for > 0, where, by employing a 
perturbation analysis, they were found to be late-time attractors. 

The endpoints of the line Li correspond to n = ±1. These represent spatially fiat 
and isotropic solutions known as the 'linear dilaton-vacuum' solutions . They are 
described by 

a = a*, 

$ = ±y2At, 

P = P*, 
(J = cr^K, 

k = 0, (A.4) 

where {a.,^., P.,a.} are constants. The "+" solution is represented by the equilib- 
rium point C in Section 3.2. 
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We now consider the cosmologies where A = in the action ( p.lO|) . A spatially 
flat solution found previously in paper II is given by 



k 



4 
- In 

0, 



M I 



3 Am ^2 



16 



16 



Am^' 



(A.5) 



where {a^,, a,,,} are arbitrary constants and time is defined over the interval 

t < (0 > 0) for the equilibrium point Si and t > (0 < 0) for the equilibrium 
point Si in Section 4.1. 

There also exists a spatially curved, isotropic solution with a constant axion field. 
It is given by 



a = 


1 




$ = 


-In 


^AMt' 


/9 = 






a = 


0"*, 




k = 


^ A 2 

-^Amo, 



(A.6) 



where {a^,, iS^^^a^,} are integration constants. The solution for t < is represented by 
the equilibrium point in Section 4.2. 

For negative Am, there are also the solutions 



a 
k 



- In 

A, 
0, 



-2AMt^ 



(A.7) 



where {a^,, a*} are constants. The ± sign corresponds to the the sign of t and 

"+" solution corresponds to the repelling equilibrium point R in Section 4.4, whereas 
the "— " solution corresponds to the attracting equilibrium point A in that Section. 
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